Quasi-condensation and coherence in the quasi-two-dimensional trapped Bose gas 
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We simulate a trapped quasi-two-dimensional Bose gas using a classical field method. To inter- 
pret our results we identify the uniform Berezinskii-Kosterlitz-Thouless (BKT) temperature Tbkt 
as where the system phase space density satisfies a critical value. We observe that density fluctu- 
ations are suppressed in the system well above Tbkt when a quasi-condensate forms as the flrst 
occurrence of degeneracy. At lower temperatures, but still above Tbkt, we observe the develop- 
ment of appreciable coherence as a prominent flnite-size effect, which manifests as bimodality in the 
momentum distribution of the system. At Tbkt algebraic decay of off-diagonal correlations occurs 
near the trap center with an exponent of 0.25, as expected for the uniform system. Our results 
characterize the low temperature phase diagram for a trapped quasi-2D Bose gas and are consistent 
with observations made in recent experiments. 
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I. INTRODUCTION 

The physics of two-dimensional systems is very dif- 
ferent from what we observe in the three-dimensional 
world. The Mermin-Wagner-Hohenberg theorem [1] |2] 
states that thermal fluctuations destroy the long-range 
order, characteristic of most phase transitions, in systems 
of reduced dimension. However, in systems that support 
topological defects, such as vortices, the existence of a 
quasi-long-range ordered (quasi-coherent) state was pre- 
dicted to occur by Berezinskii, Kosterlitz and Thouless 
(BKT) |3]|1]. The BKT superfluid transition has been ex- 
perimentally observed in liquid helium thin films [5], su- 
perconducting Josephson-j unction arrays [6J and in spin- 
polarized atomic hydrogen [7]. More recently evidence 
for the BKT transition in a dilute gas of *^Rb atoms was 
reported by the ENS group [S] ID |10]. However, strong 
fluctuations, the interplay of harmonic confinement and 
interactions, and finite size effects have made predictions 
for the low temperature phase diagram of this system the 
subject of much debate mi ttH US] [H H US]. Mean- 
field methods are inapplicable and reliable predictions 
have only recently become available from classical field 
and Quantum Monte Carlo methods |17l [TS] . 

Limited theoretical understanding has meant that the 
relationships between bimodality in the density distri- 
bution [To], algebraic decay of phase coherence [9], and 
thermal activation of phase defects |8l |9] has been un- 
clear and the subject of speculation. In particular, a 
crucial issue currently under debate is whether the criti- 
cal point identified using bimodality is distinct from 
the BKT cross-over found in |9]- In contrast to the 
ENS results a recent experiment by the NIST group 
|19) suggests that bimodality occurs at higher temper- 
ature (i.e. at lower phase space density) than the on- 
set of the BKT superfluid. In this paper we perform a 
detailed analysis of the low temperature properties of a 
quasi-2D system to explore this issue. Our results sup- 



port the NIST observations and show that bimodality oc- 
curs at temperatures above the uniform condition for the 
BKT transition |2QI |2T| . We also show measurable signa- 
tures of the degenerate components of the system: quasi- 
condensate (suppressed density fluctuations), bimodality 
(coherence), and superfluidity (algebraic decay of corre- 
lations). 
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Figure 1: (color online) Instantaneous planar density of a 
quasi-2D classical field at a range of temperatures. Inner 
white (outer black) circle marks the 1 /e-boundary of the con- 
densate (quasi-condensate) density. Vortices (+) and anti- 
vortices (-) are shown. Results indicated for various values 



of reduced temperature, T' (see Sec. II D I. Other simulation 
parameters given in the text. 
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To provide insight into the interplay of these degener- 
ate regimes, Fig. [l] shows four microstates of the quasi- 
2D trapped gas for a range of temperatures that indicate 
important features of this system: In all cases a quasi- 
condensate (as defined below) is present, and is the only 
degenerate component for the highest temperature re- 
sult Fig. [ija). As the temperature decreases [Fig. [ijb)] 
spatial coherence emerges, marked by a finite condensate 
density. However, in this temperature regime vortices are 
prolific and we frequently observe free vortices penetrat- 
ing the condensate. For the lower temperature results 
(i.e. T < Tbkt) [see Figs, [ijc) and (d)] we see that in 
the central (condensate) region density fluctuations are 
significantly reduced and vortices and anti-vortices are 
mostly found paired (i.e. in close proximity). 

II. FORMALISM 

Here we consider a harmonically trapped Bose gas de- 
scribed by the Hamiltonian 

H = j dPxi>^ i^H,p+'^^^i>^i>^i>, (1) 

where 

Hsp^^ + ^H^lx'^ + ^y + ^lz^), (2) 

is the single particle Hamiltonian, m is the atomic mass, 
and a is the s-wave scattering length. For lo^ ^ ^jJxi^y 
the z-axis degree of freedom is frozen out at sufficiently 
cold temperatures, and the system becomes quasi-2D. 
The dimensionless 2D coupling constant is g = \/Sna/az, 
with ttz — ^hjmujz- We will assume that ^ a so that 
the scattering is approximately three-dimensional [12], 
a condition well-satisfied in the ENS and NIST experi- 
ments |8]|9][10l[ig. 

In what follows we will take r = (x, y) to be the 2D 
position vector, all densities to be areal (i.e. integrated 
along z), and the trap to be radially symmetric with lo = 
ujx — LUy the radial trap frequency. Our results here are 
for the case of a quasi-2D system of *''Rb atoms with 
trap frequencies lu = 2tt x 9.4 Hz, uj^ — ^tt x 1.88 kHz 
and a coupling constant oi g = 0.107. For reference, 
in ENS experiments g « 0.13 |9], whereas in the NIST 
experiments g «0.02 [IS.. 

Our simulation method (see Appendix |A} is based 
upon a classical field (denoted by $) representation of 
the low energy modes of the system. The remaining high 
energy modes, that we refer to as the above region, are 
described by a quantum field, tp [22]. Averages of the full 
field, ^ = ^ + ip, are obtained by time-averaging the dy- 
namics of the classical field and ensemble averaging the 
properties of ■0 using a Hartree-Fock meanfield analysis 
|17l [231124] . Because the classical field approach captures 
the non-perturbative dynamics of the low energy modes, 
it is valid in the critical regime which extends over a large 



temperature region in the 2D system and was the basis 
of the definitive analysis of the uniform 2D Bose gas by 
Prokof'ev et al., in Ref. [20]. Here our approach provides 
a good description of the trapped system when: (i) g is 
small compared to unity; (ii) all modes of the classical 
field are appreciably occupied; and (iii) the fluctuation 
region is fully contained within the classical field (see 
Appendix |A]) . Our approach is computationally efficient 
and suitable for studying non-equilibrium dynamics, such 
as the vortex dynamics pictured in Fig. [l] 

We now discuss the various ways of characterizing the 
low temperature components of the quasi-2D system that 
are useful for interpreting the results of our calculations. 

A. Superfluid model of Holzmann et al. 

Recently Holzmann et al. [18] have extended a Quan- 
tum Monte Carlo method to study this system and have 
proposed a model for the superfluid component based on 
a local density application of the uniform results |20l 121) . 
In detail, they suggested a superfluid density of the form 

t 0, r > rcr 

where Vcr is the radius at which the trapped system den- 
sity equals the (temperature dependent) critical value 

n,, = A-2log(^^^ , (4) 

with C = 380 ± 3 and A = h/y/2Trmk^ {ncr>? is the 
critical phase-space density for the uniform 2D Bose gas 
to undergo the BKT transition [20]). To obtain the su- 
perfluid density using (|3j hence requires a comprehensive 
calculation of the system density (e.g. using classical fleld 
theory and Quantum Monte Carlo) from which 
be determined. 

In [ig Eq. (|3) was shown to provide a reasonable 
description of the moment of inertia obtained directly 
from Monte Carlo calculations. However this model has 
curious behavior at the transition whereby the superfluid 
density discontinuously emerges as a sharp spike when 
the central (peak) density of the system exceeds ricr- The 
temperature where this occurs is denoted TgxT |18l [25] 

US]. 

B. Quasi-condensate 

Quasi-condensate is the component of the system with 
suppressed density fluctuations, deflned by [20] 

n,,{r) = v/2(n(r))2-(n(r)2), (5) 

where n(r) = ^'^'(r)^'(r) is the density operator for the 
system. Note, that for a normal system with Gaussian 
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Figure 2: (color online) Component densities along the x- 
axis at (a) T' = 0.91 and (b) T' = 1.0 : Total density n(x) 
(dotted-line), quasi-condensate nqc{x) (dashed), condensate 
nc{x) (solid) and superfluid nsf{x) (shaded region), (c) Frac- 
tion of atoms in quasi-condensate (dots), condensate (squares) 
and superfluid (triangles). System (a): T = 30.6 nK and 
N = 9.9x 10^ system (b): T = 48.5nK and = 25.9 x 10^ 
Other parameters as discussed in the text. 

density fluctuations (ri(r)^) = 2(n(r))^ so that Ugc — 0. 
Definition (|5j was shown to furnish a universal quantity 
in the uniform gas |20] . 

C. Coherence 

There are many ways to define coherence, especially so 
for an inhomogeneous system, and here we choose to use 
a measure based on the Penrose-Onsager criterion. Nor- 
mally this criterion is that a single a macroscopic eigen- 
value of the one-body density matrix 

G(r,r') = (#nr)*(r')>, (6) 

is equivalent to the condensate number, with respective 
eigenfunction characterizing the condensate mode. For- 
mally, this analysis is only valid in the thermodynamic 
limit, and finite size effects make the unique identifica- 
tion of condensate rather fraught - especially so in the 
2D case. Here we choose to use the largest eigenvalue 
and eigenvector as a measure of coherence, and refer to 
them as condensate only for clarity |22l [27| . We denote 
the condensate density as nc(r), which provides a spa- 
tial characterization of the system coherence (e.g. the 
condensate mode boundary shown in Fig. [ij. 

D. Reduced temperature 

The classical field approach we employ uses a fixed 
classical region with its total energy as the macroscopic 
control parameter. As we vary this parameter both the 



temperature and number of particles vary (see Appendix 
[A| ) . To remove leading order effects of varying iV we scale 
all temperatures according to T' = T/Tq, where Tq is the 
quasi-2D ideal gas BEC temperature. We numerically 
calculate Tq by determining the saturated thermal cloud 
occupation according to 

N{T) = J2 (e('^-'^°)/'=^^ " ^) ' 

where {ej} are the 3D harmonic oscillator energies for our 
trap. Inverting this relation numerically we obtain To{N) 
- i.e. the temperature at which the thermal cloud of N 
atoms is saturated. The inclusion of quasi-2D effects, 
i.e. excited z-states, is typically important, producing 
an appreciably lower critical temperature than the pure 
2D prediction T2D — VONuJ^^^h/irkB |28] . 

III. RESULTS 
A. Components of the quasi-2D trapped gas 

Time-averaged density profiles are shown in Figs. [2|a) 
and|2]^b) for a system with g = 0.107. Figure |2ja) shows 
the system density when both quasi-condensate and con- 
densate components are present, whereas Fig.|2]^b) shows 
the density at a higher temperature where only a quasi- 
condensate component exists. Figure ^c) plots the de- 
generate components of the system over a large temper- 
ature range. At the highest temperatures we observe 
a small quasi-condensate fraction, which increases with 
decreasing temperature. Below T' « 0.98 the conden- 
sate fraction becomes appreciable (indicated as the co- 
herence temperature Tcoh in Fig. [2]^c)). In this regime 
we commonly see free vortices in the condensate [e.g. see 
Fig. [ijb)] which have a profound effect on its coherence 
properties. Finally, at Tbkt {T' ~ 0.93) the peak den- 
sity at the trap center satisfies Eq. (jij and the model, 
Eq. ([3}, predicts a non-zero superfluid fraction. Here 
we observe the occurrence of vortices in the condensate 
to be much less frequent than at higher temperatures, 
and those found are usually paired [e.g. see Fig. [l|c-d)]. 
As the temperature decreases further, the relative differ- 
ence between quasi-condensate, condensate and super- 
fluid fractions decreases, though they are always clearly 
distinguishable. 

An important result of this paper is that the emer- 
gence of coherence occurs at higher temperatures than 
that where the conditions for uniform system superflu- 
idity are satisfied, Tbkt, i-e. the emergence of spatial 
coherence in the trapped system occurs before the peak 
phase space density satisfies Eq. ^ . We have carried out 
simulations for g values in the range of 0.02-0.107 and 
observed similar results. 

In Fig. [3] we show the position and momentum density 
distributions for our system as a function of tempera- 
ture. In Fig. l3la) the position dependent phase-space 
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Figure 3: Density distributions versus temperature, (a) Po- 
sition phase-space density and (b) momentum space density 
distributions. The black (white) curves mark the densities at 
Tbkt {Tcoh) [also see Fig. [2]. Parameters as in Fig. |2] 



density is seen to smoothly peak as the temperature de- 
creases [29]. We have confirmed that at low temperatures 
the central density is well-described by a Thomas-Fermi- 
type profile of the form n(r) ~ n(0) — m^(jj^r^/2/i^g (also 
see [TBI). Additionally, we observe that whenever ap- 
preciable quasi-condensate is present the density profile 
bulges slightly (e.g. see Fig. |2|b)). Figure shows 
that for temperatures below T^^« 0.98 (where coherence 
is observed) that a strong bimodality is apparent in the 
momentum space density of the system associated with 
the extended spatial coherence arising from the conden- 
sate. (Note that the momentum density is shown on a 
logarithmic scale so that the change from the fiat profile 
at T' « 1 is quite profound). We emphasize that this bi- 
modality is clearly apparent at temperatures well-above 
Tbkt- 



B. Decay of coherence 

An important prediction of the BKT theory for 
the uniform system is that in the superfiuid regime 
off-diagonal correlations (phase coherence) decay alge- 
braically, i.e. G{r) ~ r~", where r is the spatial separa- 
tion and the exponent relates to the superfiuid density as 
a — l/nsf}?- At Tbkt the superfiuid transition occurs 
with UsfX^ = 4 (i.e. a = 0.25) [i] and at temperatures 
above Tbkt the system is normal with exponentially de- 
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Figure 4: (color online) (a) g^(r) for various temperatures 
(solid lines). Fitting region (light colored segments of curves) 
and algebraic fits (dashed) are also shown, (b) Peak superfiuid 
density determined from fits to algebraic decay (triangles) and 
from model density, Eq. (jsj (circles). Parameters as in Fig. [5] 



caying correlations. 

Here we investigate the nature of off-diagonal correla- 
tions in the trapped system [30]. Because of spatial in- 
homogeneity the first order correlation function depends 
on centre-of-mass and relative coordinates, and we there- 
fore choose to examine these correlations symmetrically 
about the trap center. The normalized correlation func- 
tion we evaluate is [31] 



(j>*(r/2)cI>(-r/2)) 
V(n(r/2))(n(-r/2))^ 



(8) 



and typical results for a range of temperatures are shown 
in Fig. |4]ja). We least-squares fit a model decay curve 
of the form g^{r) (x l/r" over the spatial range 1.2A < 
r < 5A to determine the exponent a. The lower spatial 
limit excludes the short range contribution of normal- 
component atoms, while the upper limit restricts the ef- 
fects of spatial inhomogeneity (typical size of cloud is of 
order 30 A). The model fits in Fig. |4|a) are poor in the 
three highest temperature cases, suggesting that the alge- 
braic fit is inappropriate. Using the fitted values of a we 
infer the superfiuid density at the trap center according 
to ns/(0) = l/a\^. The results, shown in Fig.Elb), com- 
pare well with the peak of the model superfiuid density. 
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Eq. for temperatures below Tbkt [T' « 0.93). At 
higher temperatures (where n(0) < ricr) the two results 
disagree, however in this temperature range the algebraic 
fit is poor and an exponential fit to g^{r) appears to be 
more appropriate. 

These results can provide us with an indication that 
our identification of the superfiuid transition is consis- 
tent: Using a = 0.25 as the condition for identifying the 
BKT transition and comparing to our numerical fits for 
a would suggest that Tbkt = (0.89 — 0.94)To, consistent 
with Tbkt as identified according to condition (pi. 



C. Relation to experiments 

Our results can be compared against the recent experi- 
ments at ENS and NIST. Kriiger et al. [10] (ENS) deter- 
mined the "exact critical point" by observing bimodal- 
ity in the density distribution after expansion. Those 
results suggested that the occurrence of bimodality co- 
incided with the emergence of superfluidity in the sys- 
tem. However, Clade et al. |19j (NIST), instead saw bi- 
modality occurring prior to the superfiuid transition. It 
is difficult to directly compare each experiment as they 
differ in several aspects: the ENS experiment typically 
involved several 2D systems and a large g value, whereas 
the NIST experiment isolated a single 2D system with 
a smaller g value. Our results are for a g value inter- 
mediate to both experiments, but suggest that bimodal- 
ity is distinct from superfiuidity. However, as the rela- 
tive temperature separation between these phenomena is 
quite small they may be difficult to distinguish in prac- 
tice if temperature resolution is poor. We also note that 
the ENS experimental procedure was to release the gas 
for 22 ms of time-of-fiight before absorption imaging of 
the system density. This time scale is on the order of 
the characteristic time for expansion in the weak trap 
direction {ujx « 27r x 9.4 Hz). Thus the measured den- 
sity distribution along x is in an intermediate regime not 
well-represented by the in situ spatial or momentum den- 
sities, but is in some sense a convolution of the features 
of both. A semiclassical analysis of this expansion is pre- 
sented in |3l]. While the quasi-condensate exhibits a cer- 
tain amount of spatial bulging, our conjecture is that the 
sharp bimodal feature seen in the experiment (Fig. 2 of 
jlOj l is almost certainly due to the onset of coherence [see 
Fig. |3|b)]. This could be verified experimentally using 
longer expansion times to more clearly reveal the system 
momentum distribution. 

In an earlier experiment |9] the BKT transition was 
identified by using a heterodyning scheme to determine 
when off-diagonal correlations in the central region of the 
system had an algebraic decay coefficient of a « 0.25. 
Our results in Fig. |4]show that a — 0.25 occurs where 
the peak density satisfies Eq. ([3}, which we have used 
to define Tbkt- The NIST group has developed a more 
precise procedure for measuring spatial coherence |19) . 
which should allow a more systematic investigation of 



correlations and their relationship to superfiuidity and 
bimodality. 

Another issue of concern has been the difference in 
temperature between the ideal gas condensation temper- 
ature. To, and the emergence of the degeneracy in the 
interacting system. In Kriiger et al. [10] the onset of bi- 
modality was found to occur for atom numbers "... about 
5 times higher than predicted by the semiclassical theory 
of Bose-Einstein condensation (BEC) in the ideal gas". 
This roughly corresponds to the degeneracy temperature 
being half that of the ideal gas prediction. In a subse- 
quent work the method for determining temperature was 
compared with theoretical calculations and was found to 
underestimate the true temperature by a factor of 0.6-0.7 
|34] . Although residual effects of the optical lattice po- 
tential used in the ENS experiment make computing Tq 
difficult it seems that the temperature at which bimodal- 
ity is observed is close to Tq. Recent theoretical work 
|25] has also suggested that Tq and Tbkt that are quite 
similar for the regime of the ENS experiment. Clearly 
additional work examining the critical temperature over 
a wide parameter regime will be necessary to fully un- 
derstand this relationship better. 



D. Relation to theory of Holzmann and coworkers 

Recently Holzmann and coworkers have performed 
Quantum Monte-Carlo calculations of the quasi-2D 
trapped gas jl8]. Although our results are for a system 
in a different parameter regime we find qualitatively the 
same behavior for physical predictions, e.g. the peak in 
central density curvature occurring in the transition re- 
gion (see Fig. 3 of [18]) and the appearance and behavior 
of condensate (see Fig. 4 of [H]). The Quantum Monte 
Carlo results did not analyze the transition region with 
the resolution we have in this work, and it is possible 
that the appearance of coherence before superfiuidity is 
also revealed in their results |32) . 

Holzmann and coworkers have also developed a rather 
elegant meanfield theory for predicting superfiuid tran- 
sition [25]. While providing a simple means for finding 
an approximate value for Tbkt, several aspects of this 
theory have been scrutinized in |26] and shown to be in- 
consistent. In particular, the use of (i) the total areal 
density and (ii) a temperature renormalized interaction 
strength in Eq. ^ to identify the superfiuid transition 
instead of the ground transverse mode density and the 
bare (ground transverse mode) interaction strength, re- 
spectively. Aspects (i) and (ii) were also used to identify 
the superfiuid transition in the Quantum Monte Carlo 
results and the aforementioned consistency issues may 
have contributed additional uncertainty in the transition 
temperature [18]. Because the system we consider here 
is more two-dimensional than that in [18] [33], the dis- 
tinction between total and ground transverse mode areal 
densities is negligible for our results. However, in general 
the ground transverse mode areal density is the relevant 
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quantity to be used in expression ([i} (see [26] for addi- 
tional discussion). Additionally, we note that expression 
Q is valid in the thermodynamic limit and it is not clear 
how large finite-size shifts might be. Clearly there is a 
need for more theoretical calculations in the transition re- 
gion to better clarify the relationship between coherence 
and superfiuidity. 



IV. CONCLUSIONS 

We have provided a comprehensive analysis of the 
low-temperature properties of the trapped 2D Bose gas 
and have elucidated the roles of condensation and quasi- 
condensation. A model for the superfiuid density, and 
the temperature Tbkt, was proposed in p3], and shown 
to provide a reasonable description of the system's mo- 
ment of inertia. Our results show that in the finite system 
coherence can emerge at temperatures above Tbkt- In 
this regime it is not clear if the system has some finite- 
size induced superfiuidity arising from the condensate, as 
the central coherent region is often penetrated by vortices 
[e.g. see Fig. [ijb)] which may act to destroy its superfiu- 
idity. These observations appear to be in agreement with 
the results presented in p^, and suggest that the strong 
fiuctuations seen experimentally in the initial bimodal 
state (for T > Tbxt) are associated with vortices. Our 
calculations of the off-diagonal correlations reveal that 
extended coherence emerges at temperature well-above 
Tbkt, coinciding with the appearance of bimodality in 
the momentum distribution, whereas the field develops 
algebraically decaying correlations, accurately described 
by the uniform theory of the superfiuid system, at lower 
temperatures T < Tbkt- 

As shown in Fig. [l] the classical field microstates de- 
scribe the creation of vortices and their ensuing dynam- 
ics. Future work will clarify the nature of vortex pairing, 
their typical production rates and lifetime, to better un- 
derstand their role in the equilibrium state. 
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Appendix A: REVIEW OF PGPE METHOD AND 
VALIDITY CONDITIONS 



In this work we use the classical field method known as 
the projected Gross-Pitaevskii equation (PGPE) formal- 
ism, as summarized in Refs. pTl l#l [23l [M] |30l l35l |36l l37] . 
In particular, in Refs. |171 123] this formalism is devel- 
oped for application to the quasi-2D trapped Bose gas. 



We briefiy summarise this method here before discussing 
details and validity conditions of our calculations. 



1. Use of PGPE for determining equilibrium 
properties 

A randomized classical field state ($) of definite en- 
ergy, given by the functional 



Em 



(Al) 



is constructed according to the procedure discussed in 
Sec. V of Ref . j38] . This state is then evolved according 
to the PGPE 



ih— = Hsp<^ + V 



dt 



(A2) 



where the projector, V, limits the occupation of the field 
to single particle harmonic oscillator modes below a spec- 
ified (single particle) energy cutoff ecut- The PGPE (A2l 
has been demonstrated to be ergodic so that, after the 
initial field thermalizes, time averaging can be used to ob- 
tain equilibrium quantities such as density profiles [22], 
correlation functions |22l [M] |30] . and the temperature 
and chemical potential |4H 142). 



2. Validity conditions 

The classical field method provides an accurate de- 
scription of the low energy component of the quasi-2D 
trapped Bose gas as long as three conditions can be sat- 
isfied (see sections 2 and 3 of [35] for further details) : 

(i) That g is small compared to unity: This condition 
ensures that classical fiuctuations dominate quan- 
tum fiuctuations. This condition is well-satisfied in 
our calculations {g — 0.107) and in current experi- 
ments (for the ENS group g « 0.13 [9j and for the 
NIST group g « 0.02 [19]). 

(ii) Appreciably occupied modes: In order to apply the 
classical field description we require that modes of 
the classical region are appreciably occupied, i.e. 
that the least occupied mode of the classical region 
has a mean occupation (n^i^) of order 1 or greater. 
In practice we analyze this validity requirement in 
our simulations by calculating nminj taken as the 
average the occupation of the highest energy single 
particle state (i.e. states with single particle energy 
e-fcut) M- 

(iii) Good basis: The energy cutoff ecut has to be suf- 
ficiently large that the single particle eigenstates 
provide a good basis for describing the interacting 
classical region modes. The condition for this re- 
quirement is that e^ut = ^cut ~^o^ h?gn/m, where 
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Figure 6: Relevant quantities for our classical field simula- 
tions of the quasi-2D trapped Bose gas as a function of tem- 
perature, (a) Total particle number of our results and (b) the 
respective quasi-2D ideal condensation temperature. 



Figure 5: Quantities determined from classical field evolution 
as a function of the classical field energy per particle, (a) 
Simulation temperature, (b) Average occupation of the least 
occupied mode, (c) Ratio of interaction energy scale to the 
energy cutoff. 



eo is the single particle ground state energy, and 
n is the average 2D spatial density of the classical 
field. This requirement also ensures that the fluc- 
tuation region (universal long wavelength modes) 
is contained in the classical region. This condition 
is also discussed in Ref. 1201. 



Above cutoff atoms 



Having ensured a valid classical region simulation, we 
can add the remaining low-density high-energy parti- 
cles using a Hartree-Fock description, as discussed in 
Ref. [21] . This allows us to calculate the total number of 
atoms and the total density profile for the full quasi-2D 
Bose gas. As long as the third condition above is well- 
satisfied, then these high energy atoms are well-described 
using meanfield theory [20]. 



4. Parameters for calculation 

In our calculations we have evolved classical field states 
with a cut off of Ccut = IZS^w, populated with 3750 (clas- 
sical region) ®^Rb atoms [39] for 8.5 seconds (i.e. over 80 
radial trap periods). Allowing the first 40 trap periods 
for rethermalization, we ergodically average using 2000 
sampled microstates taken over the last 40 trap periods 
of evolution. 

In Fig. |5] we show various quantities used to establish 
the validity of our calculations. We choose to keep the 
classical region size (i.e. ecut) and number of classical re- 
gion atoms fixed, and vary phase space density of the 
system by a ltering the classical region energy as given 
by Eq. (All. This change in energy directly affects the 



temperature, which we determine using Rugh's dynami- 
cal definition of temperature (e.g. see |41l H2| l. with the 
results for our parameters shown in Fig. [sj^a). 

Of the three validity conditions discussed above, only 
condition (i) can be assured a priori. Conditions (ii) and 
(iii) need to be verified a posterori as we now discuss. 
In Fig. the average occupation of the least occupied 
mode is shown. For all calculations this is larger than 
unity (and near the transition it is typically ^3 — 4) so 
that the classical field approach is well justified for these 
results. In Fig. [sjc) we show that the ratio h^gn/me'^^^ 
is small, ensuring that our cutoff is large enough to pro- 
vide a good description of the spectrum and ensuring 
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that the fluctuating region is well contained within the 
classical field. In summary Figs. |5jb)-(c) show that our 
method provides an accurate representation of the quasi- 
2D system we have simulated. In Figs. |6|a)-(b) we show 
the total atom number [Fig. iGla)] and associated ideal 
condensation temperature, ToTlFig. |6]^b)], which we use 
to rescale the temperature. Finally in Fig. |7] we show 
the thermal activation of the tight degree of freedom as 
a function of the reduced temperature. 



Figure 7: Relative thermal activation of z direction of the 
system as a function of reduced temperature 
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